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We calculate the microwave conductivity of a two band superconductor with s gap symmetry. 
Inelastic scattering is included approximately in a BCS model augmented by a temperature depen- 
dent quasiparticle scattering rate assumed, however, to be frequency independent. The possibility 
that the s-wave gap on one or the other of the electron or hole pockets is anisotropic is explored 
including cases with and without gap nodes on the Fermi surface. A comparison of our BCS results 
with those obtained in the Two Fluid Model (TFM) is provided as well as with the case of the 
cuprates where the gap has d-wave symmetry and with experimental results in Bai_ a; K a ;Fe2As2. 
The presently available microwave conductivity data in this material provides strong evidence for 
large anisotropies in the electron pocket s-wave gap. While a best fit favors a gap with nodes on the 
Fermi surface this disagrees with some but not all penetration depth measurements which would 
favor a node-less gap as do also thermal conductivity and nuclear magnetic resonance data. 



PACS numbers: 74.25.Nf, 74.20.Rp, 74.25. Fy 74.70.-b 

I. INTRODUCTION 

In a dirty isotropic s-wave BCS superconductor a so- 
called 'coherence peak' appears in the microwave conduc- 
tivity at some reduced temperature t — T/T c slightly be- 
low one. Here T c is the critical temperature at which the 
material becomes superconducting. Reducing the resid- 
ual impurity scattering rate pushes the peak closer to 
T = T c and reduces its amplitude. In the clean limit no 
coherence peak remains. Similar trends also character- 
ize the microwave response as the probing frequency is 
increased. While the weak coupling limit of Eliashbcrg 
theory reproduces the BCS results described above, an 
increase in inelastic scattering provides additional damp- 
ing effects which decrease the coherence peaks of BCS 
theory but, provided they are not too large, these re- 
main. 

A different behavior was observed in cupratc super- 
conductors very early on and this has been confirmed in 
later experiments. There is no coherence peak just below 
T c , rather a peak which can have an even larger ampli- 
tude is seen at small values of the reduced temperature 
much below T c . The peak is sensitive to residual scat- 
tering. For example doping with small amounts of Zn or 
Ni can greatly reduce the peak height and can also shift 
the temperature T at which it occurs. 1 One can under- 
stand these observations semiquantitatively within the 
phenomcnological two fluid model (TFM). Only the nor- 
mal fluid component, un(T), enters the real (absorptive) 
part of the conductivity, ai(T,cu). As the temperature is 
reduced towards zero more and more of the charge car- 
riers enter the condensate leaving less and less normal 
fluid and so decreases <J\(T, lo). But at the same time 
the quasiparticle scattering time r(T) increases and this 
increase can be sufficiently fast so as to more than com- 



pensate in oi(T, uS) for the drop in njv(T) resulting in 
its net increase with decreasing T. But t(T) can not in- 
crease indefinitely with T — > 0. Eventually it hits a maxi- 
mum value set by the finite, non-zero impurity scattering 
time even for the cleanest of samples. When this limit 
is reached, the microwave conductivity can no longer in- 
crease and, thus, will start to drop tracking the reduc- 
tion in riAr(T) as T — > 0. For a d-wave superconductor 
n n (T) is linear in T for small T but for isotropic s-wave 
it becomes exponentially small for T less than the gap 
amplitude A S (T). This clearly moves the temperature 
at which the drop in o~i(T, oS) is to be expected to higher 
values of T and should also make the drop more precipi- 
tous in isotropic s-wave than it is for <i-wave. 

In Sec. II we begin with a review of the application of 
the two fluid model to understand the microwave data 
observed in the cuprates due to the collapse of the in- 
clastic scattering at low temperatures. We also provide a 
summary of the good fit to data obtained within general- 
ized Eliashbcrg theory for d-wave symmetry and low en- 
ergy cutoff applied to the electron-boson spectral density. 
This cutoff is due to the opening of a spin gap and is iden- 
tified as the microscopic origin of reduced inelastic scat- 
tering. Next we proceed to show that the data can also 
be understood reasonably well within a simpler BCS ap- 
proach but with a phenomcnological temperature depen- 
dent but frequency independent scattering rate TgQ S (T) . 
These results alow us to understand better the limita- 
tions as well as the strength of the TFM approach. Hav- 
ing established the usefulness of the BCS formulation, 
the method is extended in Sec. Ill to the case of an s- 
wave superconductor including two bands (electron and 
hole pockets) each with a different gap value. The gaps 
could be isotropic with opposite sign corresponding to 
s^-symmetry or one of the gaps could have nodes on 
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the Fermi surface. The specific case of Bai_ 2; K 2 ;Fc2As2 
(FeAs-122) is treated extensively and comparison with 
experiment is made. Conclusions and a summary are 
found in Sec. IV 



II. MICROWAVE CONDUCTIVITY , 
PENETRATION DEPTH, AND SCATTERING 
RATES 

Many formulations of the optical conductivity start 
from a Kubo formula for the current-current correla- 
tion function. Some use a finite temperature Matsub- 
ara formalism with final analytic continuation to real 
frequencies done with Pade approximants. 2 Others pro- 
ceed within a real frequency axis formalism 3-5 for s or d- 



wave gap symmetry. For infinite free electron bands the 
integral over the energy can be performed analytically 
and the remaining integral done numerically. The in- 
put are the solutions of the appropriate Eliashberg equa- 
tions which follow once the electron-boson spectral den- 
sity a 2 F(u>) is specified. For an electron-phonon sys- 
tem a 2 F(uj) would describe the phonon exchange while 
for coupling dominantly to spin fluctuations as is envis- 
aged in the nearly antifcrromagnetic Fermi-liquid model 6 
(NAFLM) of the cuprates it describes the exchange of 
over-damped spin waves. We refer the reader to some of 
this vast literature and, here, we will not give details. 7 

The microwave conductivity of a superconductor is cal- 
culated from 
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a 1 (T,u) = ^Y^ 2v ^ / <^[f(T,v + u)-f(T,u)][A(k,u)A(k,v + u)+B(k,u)B(k,v + u)]. (1) 
k 

Here e is the charge on the electron, v the microwave frequency, /(T, uj) the Fermi-Dirac thermal occupation factor 
at temperature T and A{k,bj) and B(k, to) are, respectively, the usual charge carrier spectral density and the Gor'kov 
anomalous equivalent which is zero in the normal state. Finally, Vk x is the x component of the electron velocity at 
momentum k. The formula for the London penetration depth which is related to the zero frequency limit of the 
imaginary part of the optical conductivity by 



AY = ^o^ ,T2(T ' l/) ' (2) 



with c the velocity of light is determined by 
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The spectral densities A(k,u:) and B(k, uj) arc obtained 
in the standard way from the 2x2 Nambu matrix Green's 
function Qik, iuJ n ) with iuj n the imaginary Matsubara fre- 
quencies, iu> n = ii:T{2n + 1), n = 0,±1,±2, ... In an 
Eliashberg formulation of the theory of superconductiv- 
ity inelastic scattering is included for an electron-boson 
exchange mechanism. The superconducting gap function 
acquires a frequency dependence and the bare frequency 
uj is rcnormalizcd to uj(uj) — uj — with the 

charge carrier self energy. The static limit of the Eliash- 
berg theory reduces to BCS theory which deals only with 
impurity scattering which is static and no retardation is 
included in the pairing potential. 

As we are going to present results based only on BCS 
theory throughout the paper we need to discuss the for- 
malism applied here. It is based on the mixed symme- 
try model by Schiirrer et al. s and starts with a mixed 
symmetry order parameter A s d(0) = A s + AdV2 cos(2#) 



defining an s + d symmetric order parameter. Here, A s 
is the s-wave symmetric component and is the am- 
plitude of the d-wave component. 9 is the polar angle 
on the cylinder symmetric Fermi surface. We introduce, 
furthermore, the anisotropy parameter a by 

A s = oA , A d = v 7 ! - a 2 A (4) 
which ensures that 

A = ^(Al d (e)) g , (5) 

with (. . . )g the Fermi surface average. The gap Ao is, 
furthermore, assumed to have the standard BCS temper- 
ature dependence. According to Eq. (4) a = gives the 
pure d-wave symmetric case, while a = 1 corresponds to 
the isotropic s-wave case. The corresponding BCS equa- 
tions for the renormalized frequencies u(u>) and renor- 
malized gaps A sc i(uj) at one temperature T < T c are then 
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given in a real axis notation by 




u)(u>) 



(6a) 



A» = A s +ir^ s (^^t^=) (6b) 

/^H-A2>,0)/ e 

A d (uj) = A d (6c) 
A sd {uj) = A s {uj) + A d {Lu)V2cos(26). (6d) 

Here, Tg^g is the elastic quasiparticlc (QP) impurity scat- 
tering rate which is temperature and frequency indepen- 
dent. For convenience, we introduce 



(7) 



which, multiplied by 100, gives the percentage of the s- 
wave gap A s contained in A sd (9). This gap will have 
nodes on the Fermi surface as long as a < -\/2/3 (x < 
0.59) in a clean limit system. Nevertheless, because of 
Eq. (6b) the nodes on the Fermi surface can be lifted even 
for a < \J1j'i if the residual resistivity of the sample is 
large enough. This has also been discussed recently by 
Mishra et al. 9 



The complex optical conductivity cr(T, v) at temper- 
ature T and frequency v is calculated from the Kubo 
formula 



(3oj 



a(T,v) = ^-U ckv tanh I ^ ) [ J( w> v) - J(-w, v)\ ) , 



(8) 



with ftp the plasma frequency, (3 = l/(fcsT), and 



2J{u,v) 



1 - N(u, 0)N(lu + v,9)- P{uj, 9)P(uj + v,6) 1 + N*(u, 6)N(uj + v,6) + P*(u, 9)N(u + v, 6) 



E(u,0) + E(uj + v, 9) 



r 



E*(lu,6) - E(lu + i/, i 



where * indicates 

E(iu,e) 



the 



complex conjugate. Here, 
- A 2 sd (u; + i0+,8), N(uj,9) = 
P{u,6) = A sd {tu 



w 2 (w + 0+) 

iO+)/E(uj,6), and 
i0 + , 6)/E(lu, 9). The London penetration depth can then 
be calculated using Eq. (2) or, as was demonstrated by 
Modrc et al., 10 more conveniently in an imaginary axis 
representation of Eqs. (6) 

It is important to stress before going on to a discussion 
of BCS results as well as results based on the TFM which 
is favored in the analysis of data provided in experimental 
papers 1,11 that Eliashberg theory provides a good under- 
standing of both microwave conductivity and penetra- 
tion depth in terms of a d-wave symmetry gap function 
and an electron-boson spectral density which describes 
the coupling to over-damped spin fluctuations with a low 
frequency cutoff. This provides a temperature dependent 
inelastic QP scattering rate which can get small at low 
temperatures where it is limited only by the residual clas- 
tic impurity scattering. But the inelastic QP scattering 
is also unavoidably frequency dependent and tempera- 
ture and frequency dependence are related to each other. 
By contrast, in BCS theory the QP scattering rate t bcs 
is frequency independent as it is also the case for the 
TFM which gives as a result the temperature dependent 
scattering rate t^ m (T). In the simplest case (no ver- 



tex corrections) the optical scattering rate in the normal 
state is just twice the QP scattering rate. Nevertheless, 
one can model the inelastic QP scattering through a phe- 
nomenological temperature dependent Tg CS (T) but this 
cannot be exact as we will elaborate upon later. In the 
TFM, on the other hand, one assumes that the superfluid 
density at temperature T, ng(T), plus the normal fluid 
density ni\r(T) add up to the total electron density per 
unit volume in the normal state, denoted n. The London 
penetration depth 1/A|(T) = ikE|_ZLi22 where m is the 
electron mass. One can get the normal fluid density from 



2 n N {T) 



1 



47TAK0) 



1 



XT) 



(9) 



An optical scattering rate which we denote with t^ fm (T) 
can then be defined in terms of the microwave conduc- 
tivity as 



t tfm (T) 



c 2 1 1-A|(0)/A|(T) 



47rAl(0) 



(10) 



Before dealing with the ferropnictide superconductors 
it will prove useful to start with a very brief review of the 
situation in the cuprates. These are d-wave superconduc- 
tors but the usual isotropic s-wave Eliashberg equations 
can easily be generalized to include a momentum depen- 
dent superconducting gap which in two dimensions can 
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be taken to vary as Ao cos(2(9) with 9 an angle on the 
circular Fermi surface in the CuC>2 Brillouin zone. De- 
tails can be found in our previous papers 3-5 where we 
considered data for the penetration depth 12 and the mi- 
crowave conductivity 5 in optimally doped YBa2Cu3 0g.95 
(YBCO) single crystals and find that an excellent fit to 
both sets of data can be obtained with a spin fluctuation 
spectral form 6 (MMP form) 



a 2 F{uj) = I 2 



w/ujsf 



1 + (cj/wsf) 2 



(11) 



where / is the electron-spin fluctuation coupling strength 
and u>sf a characteristic spin fluctuation energy taken to 
be 30 meV. To fit the microwave data a low frequency cut- 
off of u c (T = 0) = 2.1 T c was applied on the MMP form 
of Eq. (11). In the NAFLM 6 this cutoff can be thought of 
as arising from the formation of a spin gap in the super- 
conducting state. This concept was already introduced 
by Nuss et al. 13 and Nicol and Carbotte 14 within the 
Marginal Fermi Liquid Model (MFLM) of Varma et al. 15 
to account for the gaping of the spin and charge suscep- 
tibility brought about by the condensation into Cooper 
pairs. The frequency cutoff u> c (T) was taken to decrease 
with increasing T according to a BCS mean field temper- 
ature dependence. The resulting normalized microwave 
conductivity a\(T) / a\(T c ) is displayed as a function of 
the reduced temperature t = T/T c as the solid (black) 
curve in Fig. 1 for the microwave frequency 34.8 GHz 
used in experiments. 1 The fit to the data is not shown 
here but it was very good. When the same model is ap- 
plied to the London penetration depth an equally good fit 
to the data reported by D. A. Bonn et al. 16 was obtained 
and it is shown as the solid (black) line in Fig. 2 for the 
normalized square of the penetration depth A|(0)/A|(T) 
in optimally doped YBCO single crystals as a function 
of the reduced temperature t. The model also provides a 
good fit to corresponding thermal conductivity data. 17 

In this paper we will use the solid (black) curves of 
Figs. 1 and 2 for microwave conductivity and penetration 
depth as representative of the cuprates and investigate 
whether or not a simpler formulation of microscopic the- 
ory, namely BCS theory, with a phcnomcnological tem- 
perature dependent scattering rate t^q S (T) can also pro- 
vide a good understanding of the data. 

We begin with a fit to the microwave conductivity 
data of Fig. 1 [solid (black) line]. The open (red) up- 
triangles are our fit with the the corresponding quasi- 
particle scattering rate denoted t^q S (T) shown as the 
open (red) up-trianglcs in Fig. 3. The fit is not unique 
and corresponds to a choice of least residual scattering 
at T = , consistent with the normalized data of Fig. 1. 
This choice is partially motivated by the recognized fact 
that the cuprates are known to be rather pure. Other 
fits all would have larger values of T^ig (T) aXT = T c as 
well as residual impurity scattering at T = 0. This am- 
biguity disappears if the plasma frequency f2 p of Eq. (8) 
is known. As a first check on the validity of our phe- 
nomenological t^q S (T) we can use it to calculate other 
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FIG. 1: (Color online) Normalized microwave conductivity 
(O / CTl Pc) vs reduced temperature t . The solid (black) line 
are the results from an Eliashberg calculation of Ref. 5 which 
fit well experimental results (Ref. 1) for a twin-free optimally 
doped YBCO single crystal. The (red) open up-triangles rep- 
resent our fit to the solid (black) curve using BCS d-wave with 
a temperature dependent scattering rate t^ s (T). The open 
(blue) down-triangles connected by a dashed-dotted line are 
additional BCS d-wave results obtained now with the scat- 



tering rate r T 



3 (T) derived from the two fluid model. 



Finally, the open (blue) down-triangles connected by the dot- 
ted line give the results of a BCS s-wave theory calculation 
using the scattering rate t~tf M -bcs(^)- 



properties. In Fig. 2 the (red) open up-triangles represent 
our results for the normalized inverse square of the pene- 
tration depth X 2 L (0)/X 2 L (T) vs T. We see good, although 
not perfect agreement with the solid (black) curve. This 
demonstrates that a BCS approach with phcnomcnologi- 
cal Tg (ig (T) fits to the microwave conductivity data can- 
not represent the penetration depth data which goes with 
it quite as accurately as we can with Eliashberg theory. 
Nevertheless, the fit is quite good and shows clearly that 
the simpler BCS approach used here can be applied with 
confidence to other systems such as the ferropnictides. 

Moreover, one can define a TFM scattering rate based 
on our BCS calculation without reference to the plasma 
frequency fl p . We define 



T TFM-BCS CO — 



1-A|(0)/A|(T) 
a[(T) 



(12) 



where cr[(T) is in computer units defined without the 
factor r2 2 /(47r) in front of the right hand side of Eq. (8). 
Results are shown in Fig. 3 as the open (blue) down- 
triangles. The points cross the open (red) up-trianglcs for 
the quasiparticle scattering rate t^ s (T) of BCS theory 
and show that the two scattering rates are not the same, 
they have a different temperature dependence. We might 
have expected them to differ only by a constant factor of 
two which is the relation expected to hold between op- 
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FIG. 2: (Color online) The normalized inverse square of the 
penetration depth Af, (0)/Ai (T) vs reduced temperature t. 
The solid (black) curve are results from an Eliashberg calcu- 
lation (Ref. 12) which fit well experimental results (Ref. 16) 
on an optimally doped YBCO single crystal. The open (red) 
up-triangles connected by a dashed line are the results of our 
BCS d-wave fit based on the scattering rate t^q S (T). The 
open (blue) down-triangles connected by a dashed-dotted line 
are the result of a BCS d-wave calculation using the scattering 
rate Ttf M _ B cs (T) derived from the two fluid model while the 
open (blue) down-triangles connected by a dotted line give 
the corresponding result of a BCS s-wave calculation. 

tical and QP residual scattering rates. But this doesn't 
hold here and shows the limitation of the concept of a 
TFM based scattering rate. We can even go further and 
use t tfm-bcs(^)/2 as an effective temperature depen- 
dent QP scattering rate in new BCS calculations. When 
this is done, we get the open (blue) down-triangles con- 
nected by a dashed-dotted line in Figs. 1 and 2 for the 
microwave conductivity and the penetration depth, re- 
spectively. The agreement with our model data [solid 
(black) line] is very poor. In particular, the peak in 
<7i (T)/<7i (T c ) is much higher in magnitude and occurs 
at higher values of the reduced temperature than in the 
model data. It is clear from this analysis that the scatter- 
ing rate obtained from a TFM analysis cannot be used in 
BCS calculations to achieve a quantitative understand- 
ing of the data. Nevertheless, it still has some usefulness 
in that it allows one to understand qualitatively how the 
collapse of the inelastic scattering at low temperatures 
can result in a peak in the microwave conductivity at 
intermediate values of the reduced temperature t. 

In Figs. 1 and 2 there is another set of open (blue) 
down-triangles connected by a dotted line. These were 
obtained in BCS calculations with t^ m bcs (T)/2 as the 
QP scattering rate in Eqs. (6a) and (6b) but now the 
gap is assumed to have s-wave symmetry (x = 1). We 
see that this assumption has a drastic effect on both mi- 
crowave conductivity and penetration depth when com- 



FIG. 3: (Color online) Scattering rates r _1 in inverse sec- 
onds as a function of temperature T for an optimally doped 
YBCO single crystal. The scattering rate t^q S (T) [open (red) 
up-triangles] was obtained by a BCS d-wave fit to the mi- 
crowave conductivity [solid (black) line in Fig. 1] and Born 
impurity scattering. Finally, the (black) open squares present 
the optical scattering rate t^ p (T), Eq. (13), at the microwave 
frequency v. Note, that in the normal state at T = 95 K the 
optical scattering rate is precisely twice the QP scattering rate 
t^q S (T = 95 K). The open (blue) down-triangles, finally, are 
the results of Eq. (12) as described in the text. 



pared with similar results for a d-wave symmetric gap. 
In particular, the microwave conductivity does not show 
a peak at intermediate reduced temperatures and is al- 
ready very small for t = 0.8. This behavior is traced to 
the much more rapid drop in normal fluid density with 
decreasing T in the s than in the d-wave case. The first 
is exponentially small at T — > while the other is linear 
in T in the same limit. This has important implications 
for the data of Hashimoto et al. 11 in Ba 1 _ :c K; C Fe2As2 as 
we will elaborate upon in the next section. 

Before turning to this discussion we make one final 
point about the pure isotropic s-wave case. In Fig. 4 
we show results of Eliashberg calculations with an MMP 
form for the electron-boson interaction as we employed 
to describe the inelastic scattering in YBCO but now an 
s-wave gap is used. The model parameters are ujsf = 
30 meV, the electron-spin fluctuation coupling strength 
I 2 was chosen to give an mass enhancement factor of 
A = 1.64, and the high energy cutoff ujd of the spec- 
trum was set to 200 meV. This resulted, together with 
the Coulomb pseudopotential /i* = 0.1 in a T c of 50 K 
more in line with the critical temperatures observed in 
the ferropnictides. Some elastic impurity scattering mod- 
eled by a constant value of the impurity scattering rate 
r hnp = 4-75 x 10 12 s" 1 was also included. We note that 
the peak in the solid curve is at lower temperature and 
is much broader than the usual coherence peak of BCS 
theory. It also reacts to the addition of elastic impu- 
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FIG. 4: The normalized microwave conductivity 
ai(T)/ai(T c ) at 28 GHz vs temperature T for an isotropic 
s-wave Eliashberg superconductor with the parameters noted 
in the figure and described in the text. The solid and dashed 
lines include residual impurity scattering T^L of 4.75 x 10 12 
inverse seconds. A low frequency cutoff in the o?F(u) 
spectrum of uj c /T c = 2.2 was applied to simulate the collapse 
of the inelastic scattering rate as T — > 0. 

rity scattering (dashed-dotted curve) in the opposite way 
to conventional BCS. 18 Increased impurity scattering de- 
pletes the magnitude of the peak, moves it to higher tem- 
peratures and narrows it considerably. These effects are 
the same as found in earlier work by Nicol and Carbotte 2 
based on the MFLM of Varma et al. 15 The physics under- 
lying the existence of the peak relates to scattering time 
variations and not to the classical coherence argument of 
BCS theory. What is important for the present paper is 
that the mechanism of the collapse in inelastic scattering 
is much less effective in producing peaks in the microwave 
conductivity for s-wave than it is for <i-wave gap symme- 
try. The fundamental difference that accounts for this 
observation is that, in s-wave the normal fluid density 
drops to zero exponentially and becomes essentially neg- 
ligible at low temperatures where the d-wave normal fluid 
density remains very significant. 

We make a final point. While we discussed two differ- 
ent scattering rates t tfm-bcs(^ 1 ) an< ^ t bcs(-^) one can 
define others which can be useful in different contexts. 
For example, the extended Drude model is often used 
to define a temperature and frequency dependent optical 
scattering rate t~^(T, uj) in terms of the complex optical 
conductivity, with 

r- 1 (T,o;) = ^Re[ ( 7- 1 (T ) a;)]. (13) 

Evaluation of Eq. (13) for the microwave frequency v = 
34.8 GHz gives the open (black) squares in Fig. 3 us- 
ing our BCS model results. It is clear that r _1 (T, v — 
34.8 GHz) is totally different from either t^ s (T) or 



t tfm-bcs(-^')- ^ u P ia Y a r °l e depending on the ques- 
tion asked. 7 ' 19,20 Finally, we would like to note that just 
above T c , r op 1 (T = 95 K) and Ttf M -bcs( t = 95 K ) a S rce 
and are twice the t^q S (T = 95 K). 



III. TWO BAND S± SUPERCONDUCTING 
STATE 

The newly discovered 21 layered fcrropnictide supercon- 
ductors display a complex band structure 22 with several 
electron and hole like pockets crossing the Fermi energy. 
Multi-band superconductivity is now well established 23 
in MgB2 which is widely believed to be a conventional 
electron-phonon mechanism superconductor with two 
bands 24-26 one with a large gap and the other much 
smaller. While in the ferropnictides there are more bands 
and the mechanism is not likely to be the electron-phonon 
interaction, 2 a minimum model that is often used is 
to include an electron band at the M and a hole band 
centered at the T points of the Brillouin zone with s- 
wave gaps of different magnitude with change in sign 
between the two referred to as s^-symmetry. 28 Angu- 
lar resolved photo emission spectroscopy (ARPES) in 
Bao.6Ko.4Fc2As2 by Ding et al. 29 gives values A2 ~ 
12meV and Ai ~ 6meV which were confirmed by 
Nakayama et al. 30 Somewhat smaller values ~ 9 meV and 
~ 4 meV are reported by the ARPES work of Evtushin- 
sky et al. 31 An important question which remains con- 
troversial is whether or not the s-wave gap on one of the 
Fermi surfaces can be sufficiently anisotropic to acquire 
a node or is it node-less. 32-36 Large anisotropics of the 
s-wave gap are certainly expected even in conventional 
electron-phonon metals. Indeed multiple plane wave cal- 
culations of the electron-phonon spectral density a 2 F(ui) 
in Pb and Al, 37-39 and first-principle calculations of the 
electron-phonon contribution to the phonon linewidth 
in Nb 40 found large anisotropies in this quantity and 
in the resulting superconducting gap values. Also in 
the high T c cuprates where the gap has d-wave sym- 
metry calculations within a BCS spin fluctuation model 
with over-damped magnons as in the work of Millis et 
al. 6 have produced gaps which go beyond the simplest 
d-wavc versions with many higher harmonics and even 
leads to mixtures of s and d-wave with profound effects 
on the resulting temperature dependence of the penetra- 
tion depth. 41-44 The penetration depth measurements in 
FeAs-122 of Hashimoto et al. 11 gave isotropic gaps with 
A2 ~ 6.8 meV and Ai ~ 3.3 meV, considerably smaller 
than ARPES but, nevertheless, implying an exponen- 
tial activated behavior at low temperatures. In sharp 
contrast, measurements by Martin et al. 45 found a non- 
exponential, close to T 2 law down to the lowest tem- 
perature measured (T ~ 0.2 T c ). While muon-spin res- 
onance experiments by Khasanov et al. 46 gives isotropic 
gaps with A2 ~ 9meV and Ai ~ 1.5 meV we note that 
this second gap is becoming rather small. Heat transport 
measurements by Luo et al. 47 are also consistent with no 
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FIG. 5: (Color online) (a) Two band s ± fit to the normalized 
microwave conductivity ai(T)/ai(T — 35 K) of Ref. 11 (open 
circles) at the microwave frequency v — 28 GHz vs temper- 
ature. In our model the small gap Ai on the hole socket is 
assumed to be isotropic s-wave while the large gap A2 on the 
electron pocket includes anisotropy characterized by the pa- 
rameter x — 0.49 fixed to get a best fit to the data [open (red) 
squares connected by a solid line] . The contribution of the gap 
tuiAi to the two band fit is shown by open (red) squares con- 
nected by a dotted line while the open (red) squares connected 
by a dashed line indicate the contribution of the second gap 
102 A2. The corresponding scattering rate i"g CS (T) is shown 
by open (red) squares in the top frame of Fig. 6. (b) The 
same as the top frame but now for a fit to experiment using a 
single anisotropic s-wave gap A2 indicated by the open (blue) 
chevron connected by a solid line. The anisotropy parameter 
is x = 0.49. 



gap nodes but they indicate that the anisotropic s-wave 
gap may be quite small in certain momentum directions. 
Finally, we mention the nuclear spin magnetic resonance 
data for 57 Fe by Yashima et a/. 48 which is consistent with 
s^-symmetry with full gaps. 

In Fig. 5(a) we show results of our two band BCS 
calculations based on Eqs. (6) for the normalized mi- 
crowave conductivity a\{T)/a\{T — 35 K) as a function 



of temperature T for a microwave frequency ^ — 28 GHz. 
The open (black) circles indicate experimental results 
by Hashimoto et al. 11 for FeAs-122, sample #3 with 
T c = 32.7 K. The open (red) squares connected by a 
solid line show the theoretical results which are seen to 
fit well the data even at low temperatures where the ab- 
sorption appears to be roughly linear in T. The corre- 
sponding temperature dependent inelastic scattering rate 
t bcs(^) * s indicated by open (red) squares in the top 
frame of Fig. 6. The small gap Ai(0) on the hole sur- 
face was assumed to be isotropic s-wave and equal to 
3.3 meV [Ai(0)/(fc s T c ) = 1.17] while the larger gap of 
amplitude A 2 (0) = 6.8 meV [A 2 (0)/(fc B T c ) = 2.4] was 
allowed to be anisotropic of the form A 2 (#) = A s d(6) = 
A s + AdV2 cos(26>). The contribution of the s-wave gap 
to A 2 is determined by the parameter x of Eq. (7) . The 
s ± -model is defined by the combination of the two gaps 
A s ± = W1A1 + iu 2 A 2 . Furthermore, the complex con- 
ductivity a(T,w) = wicrW(7» + w 2 a^(T,Lu). Here 
a^(T, to) and a^(T,u) are the complex conductivi- 
ties calculated using Eq. (8) for the two gapfunctions 
Ai(T, uS) and A 2 (T, u>), respectively. In doing so we are 
assuming that the intcrband transitions can be ignored 
which is expected since the two crW (T, u) originate from 
very different regions of momentum space. In agree- 
ment with Hashimoto et al. the weights Wi and w; 2 
have been chosen to be equal to 0.55 and 0.45, respec- 
tively. Finally, the anisotropy parameter x was allowed 
vary to get the best fit to the microwave conductivity 
data over the whole temperature range and was found 
to be equal to 0.49. Thus, we have a 49% s-wave con- 
tribution to the big gap A 2 which means that it has 
nodes on the Fermi surface although it is anisotropic s- 
wave. 8 The open (red) squares connected by a dotted 
line and a dashed line give the individual contribution 
of Ai and A 2 , respectively, to the microwave conduc- 
tivity. It is clear from this decomposition that it is the 
anisotropic gap with the nodes which contributes most 
at low temperatures as well as to the peak at ~ 17 K. 
However, the hump around 30 K is due mainly, but not 
exclusively, to the small gap. Note with reference to 
Fig. 6(a) that the temperature dependent scattering rate 
obtained [open (red) squares] is very reasonable and equal 
to 9.5 x 10 12 s _1 at T c and shows a residual scattering rate 
of Tqc S (T = 1.6 K) = 3.8 x 10 11 s" 1 indicating that sam- 
ple #3 of Hashimoto et al. 11 is rather clean, so we do 
not expect impurities to significantly alter the symmetry 
of the gap. The open circles in Fig. 6(a) are the optical 
scattering rates given in Table I of Hashimoto et al. 11 
obtained by a TFM analysis of their microwave conduc- 
tivity supplemented with their penetration depth data on 
the same sample. As we found for the case of the cuprates 
Tr^p M (T) is quite different from t^q S (T) and is smaller 
by a significant amount. The open (blue) chevrons are 
to be compared with the open (red) squares and give 
the scattering rate according to BCS theory needed to 
fit the microwave conductivity data in FeAs-122 with a 
single anisotropic gap A 2 . The fit obtained is shown 
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FIG. 6: (Color online) (a) The inelastic scattering rate 
T"g CS (T) in inverse seconds obtained from our fits to the mi- 
crowave conductivity data (Ref. If) presented in Fig. 5. The 
open (red) squares correspond to the fit presented in the top 
frame of this figure while the open (blue) chevron correspond 
to the fit presented in the bottom frame, (b) The normal- 
ized inverse square of the penetration depth Aai,(0)/A„ 6 (T) 
vs temperature T. The open (black) circles are the data for 
sample #3 of Ref. f f . The two curves indicated by open (red) 
squares and open (blue) chevron, respectively have been gen- 
erated using BCS theory the temperature dependent inelastic 
scattering rates indicated by the same symbols in the top 
frame of this figure. The open circles are from the data of 
Table I of Ref. f f . 

as the open (blue) chevrons connected by a solid line 
in Fig. 5(b) which nicely go through the experimental 
data (open circles). It is important to note that if the 
the same scattering rate t^q S (T) that fits experiment 
with the single anisotropic gap A2 is used for the small 
isotropic s-wave gap Ai we get the open (blue) chevrons 
with the dotted line through them. This curve is mainly 
confined to the temperature region above T = 20 K and 
shows, once again, that the symmetry of the gap has a 
determining affect on the temperature variation of the 
resulting microwave conductivity. 

Next we look at the temperature dependence of the 
normalized penetration depth A^ b (0)/A^ b (T) for the two 
gap s ± -model which fits well the microwave conductivity 
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T(K) 



FIG. 7: The normalized inverse square of the penetration 
depth Xat(0) / \at(T) vs temperature T. The open circles 
are the data for the FeAs-f22 sample #3 of Ref. ff . The 
solid line represents the BCS fit to the data using our two 
band s ± -model with isotropic s-wave smaller gap Ai and 
an anisotropic larger gap A2 with the anisotropy parame- 
ter x — 0.67 in the clean limit. The weights wi = 0.55 and 
to2 = 0.45. 



of Fig. 5(a). Theoretical results are shown as the open 
(red) squares on Fig. 6(b) and are to be compared with 
the open circles which are the data for sample #3 of 
Hashimoto et al. 11 Once a temperature dependent scat- 
tering rate t^q S (T) was fixed to fit the microwave con- 
ductivity there remained no adjustable parameter and, as 
one can see, the fit to the data is deficient in two ways. 
First of all, it is clear that the region near T = is not 
exponentially activated as the data indicate but is rather 
linear as one would expect from a gap with nodes. The 
data is certainly more consistent with a node-less s-wave 
gap. Secondly, the curve at higher temperatures falls 
quite a bit below the data. If we had considered a sin- 
gle anisotropic s-wavc gap fit to the microwave conduc- 
tivity data instead of our two s ± -model we would have 
obtained the open (blue) chevrons for the penetration 
depth which provide an even poorer over all fit. As we 
described in the previous section on the cuprates we do 
not expect even for a d-wavc gap that we can fit equally 
well the penetration depth with the same scattering rate 
as determined from the fit to the microwave conductivity, 
but for FeAs-122, the main problem has to do with the 
fact that the anisotropy parameter x = 0.49 obtained in 
the unconstrained fit, leads to a large gap A2 which has 
nodes. 

A natural question to ask next is what would hap- 
pen if instead of fitting the microwave conductivity we 
fit the penetration depth. In Fig. 7 we consider the BCS 
clean limit and vary the anisotropy parameter x to get 
a good fit to the data (open circles). For x = 0.67 we 
get the solid line which is in very good agreement with 
experiment. It corresponds to the mixture wjAi + W2A2 
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FIG. 8: (Color online) (a) The normalized microwave con- 
ductivity ai(T)/ai(T = 35 K) for the FeAs-122 sample #3 of 
Ref. 11 at v — 28 GHz as a function of temperature T. The 
open (red) squares connected by a solid line are the best fit 
obtained within a two band s ± -band BCS theory with a tem- 
perature dependent scattering rate t^ s (T) as to simulate the 
inelastic scattering. The smaller (hole band) gap is isotropic 
s-wave and the larger (electron band) gap is anisotropic s- 
wave with x = 0.67 chosen to give a reasonable fit to the 
penetration depth data. (See Fig. 7.) The open (red) squares 
connected to the dashed and dotted lines give the individual 
contributions of the two gaps, (b) The inelastic scattering 
rate r _1 (T) vs temperature. The open (red) squares give 
TgQ S (T), the inelastic scattering rate from our fit to the mi- 
crowave conductivity shown in the top frame of this figure. 
The open circles represent data for the FeAs-122 sample #3 
of Ref. 11. 



of the two gaps in our s -model with wi = 0.55 and 
u)2 = 0.45. We proceed with these new s ± -model pa- 
rameters to find a temperature dependent scattering rate 
to fit the microwave conductivity data. Our results are 
shown as the solid (red) squares connected by a solid line 
in Fig. 8(a). The over all fit is good including the region 
of the peak at ~ 17 K. The main deficiency is that below 
10 K where the theoretical curve drops sharply to zero 




t = T/T 



FIG. 9: (Color online) The normalized inverse square of the 
penetration depth A^ 6 (0)/A^i,(r) vs the reduced temperature 
t. The open circles represent the data by Hashimoto et al. 11 
while the open down-triangles and open up-triangles give the 
results reported by Martin et al. 45 for \ a b(0) equal to 180 nm 
and 600 nm, respectively. The open (red) squares correspond 
to our result of a BCS calculation in s symmetry where the 
large gap A2 has nodes on the electronic Fermi surface [see 
also Fig. 6(b)]. 



while the experiment still gives absorption. There arc 
two features of this fit we wish to emphasize. The open 
(red) squares connected by a dotted line give the separate 
contribution to the total from the small gap wiAi while 
the dashed curve with the open (red) squares is from the 
large, now node-less, anisotropic gap 102^2- In contrast 
to what was observed in Fig. 5(a) where the large gap 
dominated the low temperature behavior and the small 
gap contributed mainly just below T c now both contri- 
butions extend over all temperatures with the small gap 
contribution still important at lowest temperatures and 
displaying a broad peak at ~ 23 K. The BCS inelastic 
scattering rate obtained from this fit is given in Fig. 8(b) 
by open (red) squares. We note in comparison with the 
data of Fig. 6(a) that t^ cs (T = T c ) is now much bigger 
and approximately equal to 1.9 x 10 14 s _1 . This is much 
bigger than the scattering rate derived from a TFM fit 
to the data by Hashimoto et al. shown as the open cir- 
cles. While TgQ S (T = T c ) is to be interpreted as due to 
inelastic scattering it is rather big and indicates that this 
second fit to the microwave conductivity data while more 
compatible with the experimental penetration depth at 
low temperatures remains somewhat problematical. 

Figure 9 presents a comparison between experimental 
results of Hashimoto et al. 11 (open circles) and of Mar- 
tin et al. [open down-triangles for A Q {,(0) = 180 nm and 
open up-triangles for A Q f,(0) = 600 nm]. Our result of a 
BCS s^-symmetry calculation is indicated by the open 
(red) squares. They correspond to an anisotropic large 
gap A2 with nodes on the electronic Fermi surface and 
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have already been discussed in Fig. 6(b). At low temper- 
atures theory agrees well with the data by Martin et al. 45 
for A a h(0) = 180 run. For higher temperatures theory is 
always above this data set but stays significantly below 
the data for A a f,(0) = 600 nm. As A a t(0) plays the role 
of a fitting parameter in the analysis of Martin et al. 45 
a Xab ~ 400 nm will probably bring experiment closer to 
theory. Nevertheless, the low temperature dependence 
of the penetration depth as reported by Martin et al. 45 
is certainly more in line with the microwave conductivity 
data of Hashimoto et al. which shows a linear tempera- 
ture dependence of a x {T)/ai{T = 35 K) below T = 10 K 
for T — ► 0. Thus, at low temperatures there is still sub- 
stantial absorption in the system which is in contradic- 
tion to the exponentially activated behavior observed by 
Hashimoto et al. In a final point it is certainly impor- 
tant to note that the potassium content is quite different 
in the samples used in both experiments. Hashimoto et 
al. report for their sample #3 ~ 55% potassium while 
Martin et al. 45 report a potassium content of ~ 30 % for 
their sample B. 



IV. SUMMARY AND CONCLUSION 

We have reexamined the use of the two fluid model 
as a way of extracting from a combination of microwave 
conductivity and penetration depth data a temperature 
dependent scattering rate t^T pm (T) which is believed to 
model the inelastic scattering. It has the very desirable 
property that, when it is multiplied into the normal fluid 
density, it reproduces the microwave conductivity. How- 
ever, no test of its general validity has been provided. 
In this paper we take a different approach and use in- 
stead BCS theory to extract through a tight fit to the 
microwave conductivity a new temperature dependent in- 
elastic scattering rate t^ cs (T) . A comparison of Tg CS (T) 
with its TFM counterpart shows that they differ signifi- 
cantly both in absolute magnitude and in variation with 
T. This casts doubts on the quantitative significance of 
Trj7p M (T). Nevertheless, the TFM does indeed provide a 
very useful basis for a first understanding of the role of 
inelastic scattering in those phenomena. 

The fact that penetration depth information was not 
needed to extract our Tqq S (T) provides a first test of 
its validity. We use it in a BCS calculation with d-wave 
gap symmetry and find good semiquantitative agreement 
with the data in optimally doped YBCO, the same ma- 
terial used for the fit to the microwave conductivity. We 
take this small, yet significant discrepancy, as evidence 
that a temperature dependent but constant in frequency 
QP scattering rate does not capture all of the quanti- 
tative features of inelastic scattering. Because of the 
simplifications inherent in BCS theory a constant in fre- 
quency scattering rate is strictly required. But inelastic 
scattering intrinsically implies a frequency dependence 
to t -1 (T, lo) which is closely linked to its T dependence. 
For instance for electron-phonon coupling the uj 2 depen- 



dence of the spectral density a 2 F{u) at small w implies 
a T 3 law for the quasiparticle scattering rate at lu = 
and also a oj 3 law for T = 0. Similarly, for the cou- 
pling to over-damped spin fluctuations as is the case in 
the nearly antiferromagnetic Fermi-liquid model of the 
cuprates the corresponding laws are T 2 and uj 2 . Fur- 
thermore, the frequency dependence of t~ 1 (T,lu) implies 
by Kramers- Kronig transform a nontrivial (i.e.: nonzero) 
real part of the quasiparticle self energy. It is precisely 
these features that are incorporated in the ei-wave gener- 
alization of Eliashbcrg theory which we have used in our 
previous work. It is the results of our previous quantita- 
tive fit to microwave conductivity and penetration depth 
data in optimally doped YBCO single crystals that we 
have used here in our comparison with BCS. An impor- 
tant conclusion of all this is that Eliashbcrg theory is fully 
quantitative while BCS provides a good semiquantitative 
picture, its main limitation being due to the neglect of 
frequency dependence of the inelastic QP scattering rate. 



Having established the validity of BCS theory to pro- 
vide meaningful results in the present context we con- 
sidered next its generalization to the two band case. 
This is the minimum model required for a realistic treat- 
ment of superconductivity in the ferropnictides in which 
there can be several hole and electron pockets with dif- 
ferent values of the superconducting gap. A favored 
model is an s -model with isotropic s-wave gaps on 
each of the two bands and with opposite signs. We 
also consider the possibility that one of the gaps is 
anisotropic s-wave as in the work of Chubukov et al. 32 
and others. 35 Anisotropy is expected even from consider- 
ation of conventional superconductivity 36-40 and also the 
cuprates. 41-43 



A general conclusion of our analysis is that with 
isotropic s-wave it is difficult to get realistic values 
of the quasiparticle scattering rate t^q S (T) from mi- 
crowave conductivity data which shows very significant 
absorption at low temperatures as is seen in the data of 
Hashimoto et al. 11 for FcAs-122. This is also the case 
when the big gap is allowed to be anisotropic but with 
no nodes on the Fermi surface. On the other hand, if 
the gap is sufficiently anisotropic to have nodes it be- 
comes rather easy to get a fit to the FeAs-122 microwave 
conductivity data with realistic values of t^q S (T). But 
a node on one of the gaps also leads to a penetration 
depth which shows d-wave like low temperature power 
laws (i.e.: linear in temperature) rather than the expo- 
nential activation behavior found experimentally on the 
same FeAs-122 sample which is the characteristic of a 
gap with finite value everywhere on the Fermi surface. 
On the other hand, the microwave conductivity data of 
Hashimoto et al. is much more consistent with the more 
recent penetration depth data of Martin et al. 45 on a re- 
lated but not identical FeAs-122 sample. 
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